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ACCELERATION OF A CONDUCTING FLUID BY A TRAVELING MAGNETIC FIELD 


by Willis H. Braun 
Lewis Research Center 

SUMMARY 

An analysis is undertaken for the traveling-wave accelerator with the propellant 
having scalar conductivity. Two configurations of magnetic field are studied. The first 
employs a field perpendicular to the line of motion and varying sinusoidally. When the 
conductivity is assumed proportional to the density, the analysis shows that the propel- 
lant is ejected in dense pulses, there is an upper limit on efficiency of about one-half, 
and the acceleration to wave velocity can be performed in one magnetic field wavelength 
if a defined interaction parameter has a value of order one. Another accelerator using 
a magnetic wave such as might be produced by coils surrounding the channel has a some- 
what lower efficiency, which is due mainly to the loss of propellant to (or through) the 
boundaries. 


INTRODUCTION 


The desire for high specific impulses for the propulsion of interplanetary vehicles 
has led to research on several types of electromagnetic rockets. One of these is the 
traveling- magnetic -wave accelerator, which uses the principle of the induction motor. 

In this device a moving magnetic field passing through a conducting gas tends to drag the 
gas along with it. The force exerted on the gas depends on the velocity of the gas rel- 
ative to the magnetic field (see sketch). Only when the gas has reached wave velocity 
will the acceleration process stop; thus, the maximum impulse that can be imparted to 

the gas is determined by the wave speed. 

Several types of analysis have been made of devices 
of this kind. The studies of Klein and Brueckner (ref. 1) 
and Matthews (refs. 2 and 3) consider the coils that es- 
tablish the magnetic field to be pulsed very rapidly in 
succession so that the propellant is accelerated by being 


v rel. 



Force, v, 


wave 


B field 



pinched by each coil in turn. Treatments by Meyer (ref. 4) and Covert and Haldeman 
(ref. 5) regard the device as a pump that produces a steady flow against a pressure gra- 
dient. Fabri and Moulin (ref. 6) and Bernstein, et al. (ref. 7) have linearized the equa- 
tions of motion for small differences between wave and fluid speeds. 

This report attempts to preserve the nonlinear character of the equations of motion 
during an acceleration process that takes the fluid from nearly zero velocity to nearly 
wave velocity. Several simplifying assumptions about the properties of the fluid are 
made to make the analysis tractable. The analysis is applied first to an accelerator 
with the magnetic field normal to the axis of motion and operating at a high enough den- 
sity that the conductivity is a scalar; these are conditions characteristic of the experi- 
mental accelerator described by Heflinger and Schaffer (ref. 8). Secondly, there is a 
discussion of the accelerator with a magnetic field configuration similar to that of the 
device reported by Jones and Palmer (ref. 9), in which the magnetic wave has both longi- 
tudinal and transverse components. 


SYMBOLS 

A magnetic vector potential, Wb/m 

a initial coordinate of fluid particle, m 

B magnetic field, Wb/sq m 

B steady uniform magnetic field, Wb/sq m 

C cross section of accelerator, sq m 

c speed of light, m/sec 

d accelerator width, m 

F instantaneous force on propellant, N 

J conduction current, A 

j current density, A/sq m 

k wave number, 1/m 

p pressure, N/sq m 

R gas constant, J/(kg)(°K) 

r velocity ratio, v /v 

u w 

T temperature, °K 

T a average temperature, °K 
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N> 


t time, sec 

t time at which a density peak occurs, sec 

r 

u dummy variable 

v velocity of fluid particle, m/sec 

v w velocity of magnetic wave, m/sec 

X length of accelerator, m 

x longitudinal coordinate, m 

y transverse coordinate, m 

z dummy variable 

unit vector in z -direction 

a dimensionless initial coordinate, ka 

j3 magnetic field ratio, B /B 

s o 

2 

e perturbation parameter, RT /v , 

a W 

£ dimensionless transverse coordinate, ky 

77 efficiency, eq. (38) 

0 thrust, N 

k interaction parameter, a Q B q/p 0 <*> 

X wavelength of magnetic wave, m 

ji permeability, H/m 

£ dimensionless longitudinal coordinate, kx 

p density, kg/cu m 

a conductivity, mho/m 

r dimensionless time, o>t 

( p phase of magnetic wave, kx - wt 

a) angular frequency of magnetic wave, rad/sec 

Subscripts: 
o initial value 

w wave 

x,y 


vector components 



r, a partial differentiation with respect to t, a 
Superscripts: 

0, 1 zero and first order in perturbation process 
average 


CONFIGURATION 

As a simple, idealized model of the traveling -wave accelerator, consider a two- 
dimensional channel of infinite extent in the x-direction (fig. 1). An alternating wave of 
magnetic induction, which is assumed to be uniform across the channel, moves from left 
to right with speed v . Fluid enters the channel from the left at the density p Q and 
speed v Q . When the fluid reaches the origin the magnetic field begins to act upon it, 
either because it has been suddenly made electrically conducting or because it has been 
previously shielded from the field. During the acceleration of the fluid by the field, cur- 
rents flow in a direction normal to the plane of the motion and the field. The direction 
of these currents reverses as the fluid moves through a region in which the field reverses 
direction. The problem of closing the current loops at the sides of the accelerator is 
ignored here as well as the problem of the flow at the exit. If the two-dimensional ac- 
celerator is regarded as the limiting case of an annular accelerator as the two radii in- 
crease without limit, then the currents form closed loops about the axis and there is no 
problem of current-loop closure at the sides. 


MAGNETIC WAVE 

The equation that governs the vector potential zA of a magnetic field B in the x, y- 
plane (fig. 1) is 


-a- 


h A H 



Figure 1. - Schematic of traveling-wave accelerator. 
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v 2 A .J^ 

c 2 9t 2 


= -p.J 


For a quasi-static field of wave-like form, the differential operators occurring in the 
equation have magnitudes specified by 



The quantities A and J have the magnitudes 

A ~ BX 


J 


OV. 


W 


B 


2 

so that the three terms in the equation for A have the relative magnitudes 1, (v w /c) , 

5 

and Ojiv^X. Inasmuch as the desired magnitude of the wave velocity is v w ~ 10 

meters per second and the speed of light is c = 3x10** meters per second, the second 

fi 

term has a relative magnitude 10" and can be safely ignored. The third term is a 
magnetic Reynolds number based on the wave speed and length. It is assumed that the 
analysis is restricted to small values of this number: 


ajuv w X « 1 


( 1 ) 


It follows from these considerations that the vector potential is determined by 

V 2 A = 0 


If A is required to have a wave-like character, 

A = f(y)g(x - v w t) 


a solution that is symmetric about the x-axis is 
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( 2 ) 


A = — cosh ky cos(kx - wt) 
k 

For the case of a narrow channel with 


-« 1 
X 


it is nearly correct to write 


A = B Q k * cos(kx - wt) (3) 

Equation (2) represents the kind of magnetic field first introduced for traveling -wave 
machines by R. X. Meyer (ref. 4). 

LAGRANGIAN EQUATIONS OF MOTION 

Rather than to view this problem as the study of a velocity field, it proves to be con- 
venient to follow each fluid element separately - the so-called Langrangian viewpoint. 
Although both approaches are fully described only by a set of nonlinear partial differen- 
tial equations, with suitable approximations the Lagrangian equations of motion can be 
reduced to a set of ordinary differential equations. Averages over all the fluid ele- 
ments then yield mean values for the accelerator thrust, efficiency, etc. 

The Lagrangian equations of motion (Lamb, ref. 10) for a one-dimensional magneto- 
hydrodynamic flow are as follows: 

Continuity: 


_0 

at 



o 


(4) 


Momentum: 


p ^ + 9P = a / v - — \ 
at 2 9x V w dt ) 


(5) 


The continuity equation contains a derivative with respect to the distance a, which 


represents the original position of the fluid element, and, therefore, acts as a label to 
distinguish fluid elements throughout their motions. It is convenient to take -°° — a — 0 
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so that the left semiaxis (fig. 1) may be regarded as occupied by a uniform fluid that at 
time t = 0 is moved at constant speed v Q into the accelerating chamber x > 0. If C 
is the cross section of the accelerator, then the mass associated with fluid elements in 
the incremental distance Aa is p Q C Aa. During the acceleration this same mass is 
represented by pC Ax. Thus, taking the limit as Aa becomes very small yields 

P R= d ^. (6) 

p 3a 

which is a statement equivalent to equation (4) along with its initial condition. 

The momentum equation (5) shows that the Lorentz force leads to both an accelera- 
tion and a compression of the fluid. The velocity (v - 3x/3t) is the relative velocity be- 
tween the fluid and the wave. At the point where the fluid enters the accelerating 
chamber the relative velocity is greatest, and therefore the force on the fluid is maxi- 
mum. Later, as the fluid velocity approaches wave velocity, the force diminishes and 
the fluid coasts along with the field. 


NONDIMENSIONAL EQUATIONS OF MOTION 


Let the dimensionless independent variables be defined as 


2?ra 

X 


s ka = a 


Ci)t = T J 


(V 


and the dependent variable as 


£(<*, t) s kx(a, t) 


( 8 ) 


The introduction of equations (3) and (6) to (8) into the momentum equation (5) yields 

2 

3 2 ij p p 3 , p aB o /„ d%\ . 2 

3t 2 pv 2 Po da Po P°> V 0T / 


( 9 ) 


In a practical accelerator it is necessary that the terms representing the Lorentz force 
and the inertia force have the same order of magnitude. Thus, since the derivatives of 
£ are of order unity, it is necessary that 
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1 


( 10 ) 



Combined with the provision of small magnetic Reynolds number (1), this condition re- 
quires 



B 2 //i 


« 1 


( 11 ) 


Roughly, 


so that equation (10) becomes 


pv. 


w 


p v 

o 


PQ V Q V W 

B 2 /p 


« 1 


( 12 ) 


The condition (12) on the parameters of the accelerator requires that the flux of momen- 
tum in the fluid be very small compared to the magnetic pressure. If this condition is 
not met, the currents in the propellant distort the magnetic field, and the present anal- 
ysis is inadequate. 

The second term in equation (9) represents the effect the fluid elements have on each 
other’s motion through the exertion of pressure. The first factor, p/pv w , in the term 
is the square of the ratio of the local speed of sound to the speed of the wave; it can be 
anticipated that this will be very small. The second factor, p/p Q , is on the average 
about v Q /v w , also a small number. These factors should dominate the logarithmic 
derivative. Moreover, if the temperature throughout the gas does not vary too greatly 
from an average value T 0 , it is approximately correct to set p = pRT and to write 

a a 

for the pressure -gradient term (by using eq. (6)) 


pv 


A m -2- = 


2 p ca 

w 



H 

da 


-2 



An estimate of the dimensionless parameters that control the sizes of the various 
effects discussed previously can be obtained by taking the following typical operating con- 
ditions: 
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Propellant inflow, kg/sec 

Magnetic wavelength, m 

Cross-sectional area, sq m 

Magnetic field maximum, Wb/sq m 

Velocity of magnetic wave, m/sec 

Conductivity, mho/m 

Angular frequency of magnetic wave, rad/sec 

Initial density, kg/cu m 

Molecular weight, kg/g mole (argon) 

Average temperature, °K 

Under these conditions 


10" 6 
. . 1 
5X10" 3 
0. 03 
. 10 5 
. 10 4 
3?rxl0 5 
10" 5 
0. 04 
2000 


K = 


g o B o 

P G W 


« 1 


(13) 


p v v 
o w 

b 3 /p 


0.03 


(14) 


e = — — ~ 10" 4 (15) 

2 
v w 

The interaction parameter k (eq. (13)) has the appropriate size for effective accelera- 
tion. The condition of small Reynolds number as expressed in equation (12) is also met 
(eq. (14)), so there should not be large-scale distortions in the magnetic field, and 
equation (2) (or (3)) should suffice as a description of the field. Moreover, the pressure - 
gradient term has a very small coefficient e (eq. (15)); this term can, therefore, be 
expected to remain small throughout the motion. A final observation is that, under the 
conditions given previously, the self-collision time of the atoms is comparable to or 
greater than the time they remain in the accelerator. Consequently, any energy added 
to the gas during the acceleration process will not be transformed to thermal motion of 
the atoms, and the assumption of a more or less constant temperature T a should not 
cause great error. 

It will be noted that the interaction parameter k (eq. (13)) has been evaluated in 
terms of the conductivity or and density p Q at the entrance to the accelerator rather 
than locally as required by condition (10). It is assumed hereafter in the analysis that 
the ratio of conductivity to density is a constant throughout the flow, or 
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( 16 ) 


P P 0 

This assumption is made to simplify the analysis - specifically, to help reduce the 
partial differential equation of motion to an ordinary differential equation. The effect of 
this assumption on the calculated motion of the fluid will be estimated only after an ex- 
amination of the results. 

With the assumption of equation (16) for the conductivity, the magnetic Reynolds 
number for the forgoing example is pcrAv w ~ 0. 25. This estimate shows that a moderate 
amount of field distortion can be expected, and that distortion is more likely than would 
be surmised from equation (14) alone. 

It is convenient to make one more change in the dependent variable to obtain a final 
form for the momentum equation. Let 

<P(<x,t) s 4 - r (17) 


where cp is the phase of the magnetic wave at time r and position | of a particle. 
The introduction of equations (12), (15), and (17) into the momentum equation (9) yields 



dcp . 9 

+ k — sin^ffl = e 

dT 



(18) 


The initial conditions accompanying this equation are found by considering the fluid 
particle that arrives at the origin (£ = 0) at time t. Its particle number (initial location) 
is 

a = - v t 
o 


or 


a = - rr 


r s V v w 


(19) 


Equation (19) describes the initial curve in the a, r-plane, along which the phase ip 
and its derivatives must be given: 


(p = - T 


a 


( 20 ) 
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( 21 ) 


dcp _ 9| 
3 t 3 r 


1 = r - 1 


dcp _ d£ _ Po 
da da p 


( 22 ) 


In equation (21), 3|/3r is the fluid velocity normalized to the wave velocity. The con- 
tinuity equation (6) has been employed in equation (22). Equations (18) to (22) comprise 
the complete boundary- value problem for the system of fluid particles which is accel- 
erated by the magnetic wave. 


PERTURBATION ANALYSIS 

Inasmuch as a small dimensionless number e arises naturally in the formulation of 
the equation of motion (18), it is convenient to regard it as a perturbation parameter for 
reducing the partial differential equation to a system of ordinary equations. If the de- 
pendent variable cp is expanded in powers of e, 

cp = qp + ecp^ + . . . (23) 


the coefficients of the first two powers in the momentum equation (18) are as follows: 
.0 


For e 


d 2 (p° d(p° 2 0 

— — + k -Z- sin z cp u = 0 

3t 2 9t 


(24) 


For e : 


.2 0 

3 cp 


= 2.1 


+ * i- W 1 Sin 2 <p°) , 

* 2 ST 

\ da , 


(25) 


The corresponding initial conditions are 
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II I II HI 


( 26 ) 


(27) 


The zero-order equation of motion (24) represents the motion of one particle alone 
that is unaffected by the presence in the channel of other particles; this is the result of 
the disappearance of the pressure-gradient term to this order of magnitude. Physically, 
the meaning is that the magnetic wave moves so fast (with speed v ) and accelerates the 
fluid so quickly that the forces due to pressure gradients, which propagate with the 
speed of sound, do not have time to readjust the fluid distribution during the acceleration 
process. Of course, this condition will no longer hold after the fluid has been accel- 
erated to nearly the wave speed, for then (dcp/dr) ~ 0 and the Lorentz force term is 
negligible. The effects of pressure gradient will then take over. 

The first-order equation (25) accounts for the effect of pressure gradient; that is, 
q > * makes a correction on qp to account for the neglect of the pressure term in equa- 
tion (24). Moreover, both equations (24) and (25) can now be regarded as ordinary dif- 
ferential equations in r because the partial derivatives of any qP with respect to the 
material coordinate a are always transferred to the equation of order n + 1. The 
initial conditions (eqs. (20) and (22)) are satisfied entirely by qP. The higher order cp n 
have homogeneous initial conditions. 

The first integral of the zero-order equation is (when using the notation <p ® for 

dqP/dr) 

qp{ t, a) = r - 1 - — (l qP - 2 — - sin 2 qp + sin 2 — j (28) 

7 4 y r r J 

The variable r may be found as a function of qp by integration from the initial condi- 
tion (26). Then 


I 
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(29) 


r = 



dz 

<?z(z; «) 


4-i_ 

Equation (29) specifies the time r at which the a 1 fluid element arrives at the 
location in the wave having phase qP . The corresponding position and velocity are, to 
zero order in e, 


i = 


0 

<P + T 


(30) 


and 


v = v £ = v (qp + 1) (31) 

w^r w vr r ' ' ' 

The only remaining piece of information to be found is the density. When equation (29) is 
differentiated with respect to a and when the expression for the density (eq. (6)) is 
used, there results to lowest order in e 



This expression for the density, when evaluated on the initial curve, conforms to the 
initial condition (eq. (22)). 

The first-order equation (25) with initial conditions (eq. (27)) can also be integrated 
directly. The first integral is 


(p* + Kcpl sin 2 qp 



Vaa^’ 01 ) dz 



and the second integral is 


(33) 
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<pH<P°i a ) = exp 



sin^ z dz 

<Pz ( z ; «) 




dz 

<?z(z; «) 


( 34 ) 


For the explicit evaluation of the integrals in equations (33) and (34), it is necessary to 
have cp^ xa as well as the expressions (eqs. (28) and (32)) for cp® and cp ® . Differentia- 
tion of equations (28) and (32) yields 



To first order, the position, density, and velocity of a particle can now be written 
in terms of cp ^ and its derivatives as 


4 = 


< p® + ecp^ 


+ r 


p o o 

j=<P a + ecp 


1 

a 
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V , o 1 

— = l + cp T + e<p T 

w 


EFFICIENCY AND THRUST 

Energy is being added to the accelerated fluid in two ways: the kinetic energy is 
being increased by the acceleration process, and heat is being added by Joule dissipation. 
Since the propellant atoms are isothermal, the heat added must manifest itself in excita- 
tion and ionization of the fluid and in electron heating. When it is assumed that the heat 
added is sufficient to maintain the conductivity of the propellant and when the loss of 
kinetic energy by drag at the walls is neglected, the efficiency is 


V = 


kineti c energy added 

kinetic energy added + Joule heat added 


( 36 ) 


The kinetic energy added per unit mass is 


I (v 2 - v 2 ) =i v 2 ({ 2 - r 2 ) v 2 [(,? + l ) 2 - r 2 + *2<e\<.<e° T + « + •••] 


The Joule heat per unit mass, J /crp, when integrated over a particle trajectory is 

y 2 2 / sr ° 

cp T sin cp dcp = v k 

i /r J a /r 


c - 

''-a/r P 


- (v - v) 2 B 2 dt 

- ' w ' 


= V 2 k f (P sin 2 cp d<p = v 2 k f <p Q sin 2 <p° dcp° 
Ja/r Ja/r 


+ ev^/c f sin 2 <p° + cp^cp® sin + <p® sin 2 (p ® -^-^d<p® + ev^Kcp^cp® sin 2 cp ® 

W *4/r V 9 ^ / 


( 37 ) 


To zero order in e, 


15 


V - 


, 0 iN 2 2 

(<P r + 1) - r 


/ 0 ,\2 o /'(» 0 0-2 0,0 2 

(<p +1) + 2/c / ^ <p sin cp dcp - r 
T -Ax /r T 


After the integration in the denominator is performed, it is found that 


7 ] = 


<pz. + 1 + r £ + r 

T _ T 


( 38 ) 


The efficiency to order zero associated with a fluid particle is simply the average of the 
final and initial dimensionless velocities of the particle. Since the dimensionless veloc- 
ity has an upper limit of unity and r is very small, the efficiency has an upper limit of 
about 1/2, a result previously obtained by other investigators (ref. 8). 

An average efficiency may be obtained for any station along the accelerator by 
averaging over all particles that enter during a half cycle of the magnetic wave. Thus, 


where 





The thrust of the accelerator is calculated by finding the total force acting on the 
fluid at any one time and then averaging over a period. If X is the length of the accel- 
erator and C is its cross section, the force acting at any time is 


pc , 2-1 /'“< x )/ r . ; 

F = X ”('» -“(B C dx = - <r 0 v 0 B 0 Ck / /r= T Vr sm 



The average of this force over a period of the magnetic field is the thrust 



F(r)dr 
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Because the force is periodic it is permissible to change the area of integration in the 
O', r-plane to represent the integration over the paths of a set of particles which come 
into the accelerator during a half cycle of the magnetic field. Then 



Equation (28) has been used to obtain the last integral. If the definition of k is 
employed, then 

® = VoW - r > 

Thus, the thrust is just the mass flow into the accelerator times the average increment 
of specific impulse, a result which bears a close resemblance to that for a steady flow. 

It should be emphasized, however, that the average appearing here is not a time average, 
but an average over fluid elements of equal mass. 

NUMERICAL RESULTS 

The integrals appearing in equations (29) and (32) for the zero-order motion of 
50 fluid elements in a wave half cycle have been performed. These particles enter the 
accelerator at uniformly spaced times of -a/r = 0, 7t/50 . . . tt. Along with equa- 



0 2 4 6 8 10 12 14 16 18 20 22 24 26 

Dimensionless longitudinal distance, £, =-2tkI\ 


0 12 3 4 

Number of wavelengths 

Figured - Average efficiency of traveling-wave accelerator. Velocity ratio, v Q /v w , 0.01. 
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tion (28) these integrations make it possible to compute the zero-order velocity, effi- 
ciency, mass density, and current density associated with each element of fluid. Inas- 
much as efficiency and thrust are linearly related according to 

0 = 2p o v o v w C( ^ - r > 

it is sufficient to work with the efficiency. 

Figure 2 shows how the efficiency improves with length of run down the channel for 
various interaction parameters k. The waviness in the curves can be attributed to the 
periodic nature of the driving force. Generally, an accelerator can be expected to be an 
integral number of waves in length in order that the driving coils can properly establish 
the traveling wave. For a one-wavelength accelerator, a strong interaction parameter 
such as /c = 1.0 seems to be greater than required to successfully accelerate the fluid 
in the single wavelength. At the other extreme, a value k = 0. 1 does not provide enough 
interaction between wave and fluid to get the fluid near the wave speed at the exit. For a 
one-wavelength accelerator, it appears that some intermediate value of the interaction 
parameter, say k = 0. 3 to 0. 5, would be appropriate. 

Figure 3 shows how the field, velocity, current density, and mass density vary in 
time at various stations and interaction parameters. The most striking features of these 
distributions are the sharp accumulations of mass and current. The peaking of these 
quantities becomes more severe as the fluid moves down the accelerator. Only when the 
fluid has nearly reached the wave velocity will pressure effects, as represented in first- 
order calculations, take over and redistribute the fluid in a more uniform manner in the 
beam. The ratio of the mass density in the peak to the mass density midway between 
peaks can vary over several orders of magnitude; values range from 30 to 3000 in the 
numerical examples computed. The fluid elements, which enter the accelerator chamber 
coincident with a node in the magnetic wave, form the centers of the regions of high den- 
sity. 

Although the first-order function cp ^ has not been calculated, the right side of 
equation (25), which determines its magnitude, has been computed at every step in the 
calculations. This term, which may be regarded, alternatively, as the neglected term 
of the zero-order equation (24), is generally of order unity, except possibly when a fluid 
element enters the accelerator in which case it may be as high as 90. This is still well 
within the allowable range set by the estimate of e (eq. (15)). Therefore the zero-order 
solution should be a good approximation to the complete solution; the pressure effects 
are not really important. 

A comparison of the present work with the computations of Fabri and Moulin (ref. 6) 
is made in figures 4 and 5. There is a limit to the aptness of such a comparison because 
the analysis of Fabri and Moulin is valid only when the difference between initial fluid 
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Mass density ratio, p/p 0 Current density, Velocity ra- Field ratio, B/B 0 Mass density ratio, p/p 0 Current density, Velocity ratio, Field ratio, B/B 






Mass density ratio, p/p 0 Velocity ratio, v/v 



velocity and magnetic wave velocity is small, whereas the present analysis holds when 
the magnetic wave velocity is much greater than the speed of sound and, therefore, the 
inlet velocity. The present analysis, consequently, drops the pressure-gradient term. 
The example used for comparison is in the range of inlet velocities for which the assump- 
tions at Fabri and Moulin are valid. Nevertheless, the results of the two computations 
show much similarity. Figure 4 shows maximum and minimum velocities observed at 
ainy station along the accelerator, and figure 5 shows maximum and minimum densities. 
The two solutions yield much the same variation of the extreme values near the inlet, 
but further downstream the effect of the pressure-gradient term is very apparent in the 
solution of Fabri and Moulin. The relatively narrow bands of velocities and densities 
in that solution must be due principally to the tendency of the pressure gradient to keep 
the flow uniform. Although the time-average velocity has not been computed in the 
present computations, it appears from figure 4 that the average velocities for the two 
methods should be in good agreement. 

The curves of figures 4 and 5 also show that the maximum velocities and densities 
for the two methods differ more in the cited example than do the minimum values. This 
fact may be associated with the assumption in the present work that the conductivity is 
proportional to the density, as contrasted to the assumption of Fabri and Moulin of con- 
stant conductivity, for the former assumption leads to regions of high coupling between 
fluid and wave wherever increases in density occur. Eventually a large fraction of the 
fluid accumulates in high-density regions which are moved close to wave speed. If the 
minimum velocities and densities for the two methods are generally in more agreement, 
as is indicated by this one example, then the present method will give a higher estimate 
of the efficiency of the acceleration process than would be expected under an assumption 
of constant conductivity. 

In considering the previous example as described by the two methods of calculation, 
it should be kept in mind that the example is not a very realistic one because of the small 
differences in velocity between wave and fluid. It is chosen to make comparison easy 
with the method of Fabri and Moulin. The present method, valid for large velocity dif- 
ferences, quite correctly ignores pressure gradients under practical conditions. 

One other question that arises in connection with the assumption of conductivity pro- 
portional to density is whether the Joule heating can provide sufficient energy to keep the 
gas in a state of high conductivity as the postulate requires. A simple calculation suf- 
fices to show that it probably can. From equation (37) the energy dissipated per unit 

o 

mass over the time the gas is in the accelerator is of order v w /lO. For argon this is 
equivalent to about 40 electron volts per atom. Since most of the dissipation occurs near 

_ O 

the entrance and the decay time for free electrons in argon is 10 second as opposed to 
10 ® second acceleration time, there appears to be enough energy to keep a sufficient 
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fraction (say 1 percent) of the atoms ionized and to raise the electron temperature, even 
if excitation of atoms absorbs a large amount. 


LONGITUDINAL -WAVE ACCELERATOR 

Another form of traveling-wave accelerator uses a moving magnetic field with a 
potential 

. B o . 

A sin cp sinh ky 

k 

as contrasted to equation (2). The field resembles qualitatively that produced by coils 
wrapped around the channel at equal intervals Ax and having equal phase differences. 
In addition to the traveling magnetic wave, a uniform field B parallel to the direction 
of motion may be impressed to help contain the fluid and keep it from the walls. The 
equations of motion are, neglecting pressure gradients, 

^ - (i - » W )B* + yB x B y 

J y = < x ' v w> B x B y - 

and the initial conditions on each fluid element are 

"cp = r - 1 

1 = 0 

The components of magnetic field are 

B x = - B q sin cp cosh ky + B g 

By = B q cos y sinh ky 

In terms of the definitions, 
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£ = ky 

£, 

£ - ]3)cos cp sinh £^ 

£ = - /cj<p(sin cp cosh £- j3)cos cp sinh £ + £(sin cp cosh £ - /3) 2 
The efficiency for a particle is given by 


the dimensionless equations of motion are, for small 

•• r 2 2 

<p = -Kup cos cp sinh £ + £(sin cp cosh 


(39) 

(40) 


V = 


, . 1N 2 2 
(<P + 1) ~ r 

f + (* + l ) 2 + £ 2 - r 2 


where 


H = kvI j cos cp sinh £ + £(sin <p cosh £ - /3)J 2 dr 

The equations of motion (39) and (40) have been integrated for two cases - one with a 
confining field (/3 = 1) and one without 0 = 0). In each case the calculations have been 
terminated at a lateral position designated as the ”wall M ; that is, no particles are con- 
sidered which enter at a distance farther fio m the center line than the wall position, and 
the trajectory of any particle is no longer followed after it crosses the wall ordinate. 

In the first of these examples, the width of the accelerator is about x/10. Trajec- 
tories have been obtained for fluid elements entering at both the wave node and the wave 
crest. The lateral displacements of the elements at the entrance are £ = 0. 01, 

0. 02 . . ., 0. 29, 0. 30. For the second example, an accelerator with no confining field, 
the width is approximately 0. 6 x and the entrance ordinates are £ = 0. 1, 0. 2 . . . , 

1. 8. The width, in this case, corresponds to the optimum coil radius of 0. 3 X found by 
Palmer, Jones, and Seikel (ref. 11) for the efficient production of a traveling wave by a 
system of equally spaced circular coils. In the second example, also, trajectories for 
fluid elements entering at the wave crest and node have been found. For both cases the 
interaction parameter is k = 1 and the velocity ratio is r = 0. 01. 
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Figure 6. - Propellant trajectories in longitudinal wave accelerator. Interaction parameter, k, 1; magnetic 
field ratio, B^B^ 1; velocity ratio, v Q /v w , 0. 01. 



Dimensionless longitudinal distance, £ = 27tx/A 


0 .1 .2 .3 .4 

Number of wavelengths 

Figure 7. - Loss of propellant entering at wave node to walls of longitudinal wave 
accelerator. Phase of magnetic wave, <p, zero when dimensionless longitudinal 
distance, is zero; accelerator width, approximately one-tenth wavelength of 
magnetic wave. A; interaction parameter, k, 1; magnetic field ratio, Bg/B^ 1; 
velocity ratio, v Q /v w , 0. 01. 

Figure 6 shows propellant trajectories in the first case when there is a confining 
field. It is apparent that eventually all fluid elements will strike the walls. A large 
fraction is lost even on the first undulation, which occurs within one-tenth of a wave- 
length from the entrance. This effect is displayed better in figure 7, which shows the 
fractional loss as a function of accelerator length. It is apparent that the inclusion of a 
moderate (/3 = 1) confining field cannot prevent the loss to (or through) the boundaries of 
the accelerator of a sizable fraction of the fluid. The efficiencies associated with the 
fluid elements at the time they hit the wall are 15 percent or less. Calculations for 
fluid elements entering at the wave crest show essentially the same behavior (fig. 8), 
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Fractional loss 


except that they are thrown to the walls 
roughly a quarter of a cycle earlier. 

When the accelerator channel is wide, as 
specified in the previous second example, the 
trajectories of the fluid elements show a con- 
siderable variation, which depends on the 
entrance coordinate £. Some examples are 
given in figures 9(a) and (b). The time that 
each particle is in the accelerator before 
striking the wall can be estimated by the num- 
ber of undulations it undergoes. Each rise or 
fall represents the passage of one -fourth 
wavelength of the magnetic field and, there- 
fore, a quarter period of time. Although a 
fluid element close to the centerline may 
reach the wall at an earlier axial position 
than an element further out, it does so at a 
later time. An interesting feature of the tra- 
jectories shown is that they all terminate in a 
portion that is nearly rectilinear motion. 
Moreover, this motion is found to take place 
at nearly constant velocity. In this last portion of its flight, each fluid element affixes 
itself to a magnetic field line and rides out of the accelerator at constant velocity, much 
as was the case in the normal- field accelerator, but in this case the velocity is not paral- 
lel to the accelerator axis. 

Figures 10 and 11 show the fractional loss to the walls for fluid entering at the wave 
nodes and crests, respectively. There are large losses to the walls again; however, the 
average efficiency of the particles is higher than in the narrow channel case. For ex- 
ample, the fluid elements that reach an axial position of one wavelength (ij = 2 tt) have an 
average efficiency of about 0. 4, and, if no contribution is allowed from those elements 
which strike the wall, the average efficiency for all is about 0. 26. This is a consider 
able improvement over the previous case and reinforces the findings of Palmer, Jones, 
and Seikel (ref. 11) that the wider accelerating channel should be the more efficient one. 

Finally, it appears that a confining field is less effective in containing the propellant 
than is the widening of the channel. A comparison of figures 7 and 10 and figures 8 
and 11 shows that most of the propellant in the narrow accelerator having a confining 
field is lost within a wavelength, whereas in the wider channel having no confining field a 
correspondingly high level of loss is not reached for several wavelengths. 
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Figure 8. - Loss of propellant entering at wave crest to 
walls of longitudinal wave accelerator. Phase of magnetic 
wave, <p, 7t/ 2 when dimensionless longitudinal distance, 
F,, is zero; accelerator width, approximately one-tenth 
wavelength of magnetic wave, A; interaction parameter, 
k, 1; magnetic field ratio, B s /Bq, 1; velocity ratio, 

W 0.01. 
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(a) Particles entering at maximum longitudinal field. Phase of magnetic wave, <p, nl2 when dimensionless longitudinal distance, £,, 
is zero. 

Figure 9. - Propellant trajectories in longitudinal wave accelerator. Accelerator width, approximately six-tenths wavelength of mag- 
netic wave, A; interaction parameter, k, 1; magnetic field ratio, BJBq, 0; velocity ratio, v 0 /v w , 0.01. 
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<b) Particles entering at minimum longitudinal field. Phase of magnetic wave, c p, zero when dimensionless longitudinal distance, 
is zero. 

Figure 9. - Concluded. 
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Figure 10. - Loss to walls of propellant entering at minimum longitudinal field. Phase of 
magnetic wave, <p, zero when dimensionless longitudinal distance, £, is zero; accelerator 
width, approximately six-tenths wavelength of magnetic wave, interaction parameter, 
k, 1; magnetic field ratio, B s /Bq, 0; velocity ratio, v Q /v w , 0.01. 
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Figure 11. - Loss to walls of propellant entering at maximum longitudinal field. Phase of 
magnetic wave, <p, irl2 when dimensionless longitudinal distance, is zero; accelerator 
width, approximately six-tenths wavelength of magnetic wave, interaction parameter, 
k, 1; magnetic field ratio, B^Bq, 0; velocity ratio, v Q /v w , 0.01. 
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CONCLUDING REMARKS 


Principally, the limits on the accuracy of the analysis made here come from the 
assumption that the conductivity is proportional to the density. The assumption was 
made for convenience in reducing the equation to manageable form, and it is difficult to 
say just how the conductivity should be represented in an unsteady, nonuniform situation 
as is presented in the accelerator. The most that can be said of the effect of the con- 
ductivity assumption is that it probably exaggerates the variations of density and velocity 
that occur in the flow, although the average values should be representative. In this con- 
nection, the upper limit of 0. 5 on the efficiency is a good indication of the performance 
of the normal- field accelerator. 

The solutions to the longitudinal-wave accelerator have a further limitation in the 
amount of information they can offer because it is not possible to obtain the density with- 
out first finding the trajectories of a very great many fluid elements. Nevertheless, the 
relative efficiencies for the two calculated cases - narrow and wide channel - are suffi- 
cient to illustrate that accelerators of greater ratio of coil radius to wavelength are more 
efficient, which is just a reflection of the increased efficiency of acceleration of a fluid 
element with the starting value of y or £ , and the attendant increase in the normal 
component of field. Even the wide-channel accelerator of this type, however, attains 
only about one-half (0. 26) the efficiency reached by the normal field configuration, and 
then only if the accelerator is several wavelengths long. 
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